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Abstract 

Consider a fluid flowing through a junction between two pipes with 
different sections. Its evolution is described by the 2D or 3D Euler 
equations, whose analytical theory is far from complete and whose 
numerical treatment may be rather costly. This note compares different 
ID approaches to this phenomenon. 
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1 Introduction 

The motion of a fluid in a pipe is described by the one dimensional p-system 
in Eulerian coordinates, which reads 

= 

= 0. C-i) 



where p is the fluid mass density, q its linear momentum density, p its 
pressure, x is the space coordinate and t is time. Assume that the section 
of the pipe has a change due to a junction, sited at, say, x = 0. The motion 
of the fluid can then be described by (|1.1|) together with a condition on the 
traces of the thermodynamic variables at x = 0, i.e. a condition of the type 

$ {ai,p{t, 0-), q{t, 0-); a„ 0+), q{t, 0+)) = ; (1.2) 
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here ai and are the pipes' sections to the left and to the right of the 
junction. Various choices of the function <I> are present in the hterature, 
see for instance [21 |3l El El EJ [12]. Here, we compare the various definitions 
on the basis, a priori, of their physical meaning and, a posteriori, of their 
analytical properties. 

This phenomenon is essentially 3D, due to turbulence, or may be re- 
duced to 2D, but it is intrinsically multi-dimensional: "It is probable that, 
in the future, changes of cross-sectional area in pipework systems will be in- 
corporated as multi- dimensional patches in one- dimensional simulations. 
[201 Section 6.1., p. 262]. However, the numerical treatment of 2D or 3D 
hyperbolic systems is far more expensive than that of ID systems. We refer 
to [131 Table 1] for a striking comparison between CPU times in 2D and 
ID numerical integrations of Euler equations in the case of a "T" junction. 
Besides, as is well known, basic analytical questions, such as well posedness, 
about 2D or 3D Euler equations are still unanswered. Aiming also at the 
case of more complex gas networks, the availability of ID simple, though 
approximate, "patches" (with reference to the citation above) may be of 
great help. All this implies, in turn, to choose a specific function $ in (|1.2|) . 
Below, we compare various choices appeared in the recent literature. 

First, we deal with solutions to Riemann Problems for (|1.1|) at a junction 
between two pipes with different sections. By Riemann Problem for (jl.ip 
we mean 



Problem plays an essential role in the construction of solutions to the general 
Cauchy Problem for Analogously, the solution of ([ri])-([r2]) at 

a junction is essential for the development of the theory of the p-system on 
general networks, see |llj . 

Following various results in the literature, such as [2l[3l[5l[6l[71[8l[9l El 
[TC] . we restrict our attention to subsonic fiows in both pipes. 

Remark that the modeling of compressible fiuid flowing through pipes 
with sudden enlargements can be achieved through several entirely different 
techniques, for instance the isothermal lattice-Boltzmann methods, see [I], 
or the full system of Navier-Stokes equations, see [18], whereas experimental 
data are also available in the literature, see [10 ^ 1161 [T7] . 

The next section is devoted to ()1.3p and to the Cauchy problem for (jl.ip - 
()1.2p . both in general and for specific choices of the junction condition 
Then, Section [3] displays some numerical solutions. All the analytical details 
are gathered in Section jH 



dtp + d^q = 0, 

dtq + d^{q'^/p+p{p)) = 0, 



(1.3) 




where {pi,qi) and {pi,qi) are fixed in M"*" x M. Recall that the Riemann 
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2 Analytical Results 



Throughout, we identify u = (p, q). On the pressure law p, we assume 

(EoS) p G C2(M+;M+) is such that for all p > 0, p'{p) > and p"{p) > 0. 

The classical example is the 7-law, where p{p) = k (p/p)"' , where 7 > 1 and 
p > 0. Recall for later use the expression of the flow of the linear momentum: 

P{u) = ^+p{p). 
P 

Throughout, the present analysis is limited to the subsonic region given by 

Ao = |uGi+xM: Ai(u) < < A2(n)| , (2.1) 

where Aj is the i-th eigenvalue of (jl.ip . see ()4.ip . Recall the following defi- 
nition of solution to the Riemann Problem (|1.3p under condition (|1.2p . see 
also P Definition 2.2]. 

Definition 2.1 A weak ^-solution to the Riemann Problem is a map 

(2.2) 



u{t) G BV(R;M+x]R) for a.e. t e 
such that 

1. for all If G C^(M^ x ]R;M) whose support does not intersect x = 




P 
Q 



Q 

Pip,q) 



dxf dxdt = 0; 



2. for a.e. x G M, ti(0,x) coincides with the initial datum in hl.S^) : 

3. for a.e. t G M"*", the junction condition at the junction is met. 

We consider the following properties of the junction condition p.2p . which 
we rewrite here as 

^{ai,ui]ar,Ur) = . (2.3) 

(*0) Regularity: $ G (^(M+ x Aof;R'^y where ^0 is given by Oil) . 

Moreover, the 2x2 matrix Du^^{ai,ui;ar,Ur) is invertible, for all 
a/ , Or > and ui,Ur G Aq. 

($1) No-junction case: for all a > and ui,Ur G Aq, ^{a,ui, a,Ur) = if 
and only if ui = Ur- 
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($2) Left-right symmetry: for all ai,ar > and {pi,qi),{pr,qr) £ ^o, 
^{ai,pi,qi;ar,pr,qr) = if and only if ^{ar,Pr, -qr^a^pu -qi) = 0. 

($3) Consistency: for all positive ai,am,CLr and all ui,Um,Ur G ^o, if 
^{ai,ui;am,Um) = and ^ {am, Um;ar,Ur) = then $(a;, u;; a^, u^) = 
0. 

($4) Hydrostatic limit: for all positive ai,ar and for all densities pi,Pr, 
^{ai,pi,0;ar,Pr,0) = if and only if p{pi) =p{pr)- 

Moreover, by an immediate extension of [9, Lemma 2.1], ($[0]) ensures 
that (12. 3p implicitly defines a map = T{ui;ai,ar) in a neighborhood 
of a subsonic state satisfying ^{ai,ui; ar,Ur) = 0. In turn, this implies the 
local well posedness of Cauchy problems for data near to stationary solu- 
tions; see [3 |8l [9] and Proposition 12.21 below. Here, "near" is meant in the 
sense of the total variation. 

Proposition 2.2 Let p satisfy (EoS), <I> satisfy (^D^). Then, for all posi- 
tive ai,ar and ui,Ur S such that ^{ai,ui;ar,Ur) = 0, there exist neighbor- 
hoods Ai,Ar of ai,ar and Ui,Ur of ui,Ur such that for all ai £ Ai, Ur £ Ar 
and ui GUi, there exists a unique Ur € Ur such that the map 

is a stationary weak ^-solution to (T7Jl)-^K^ in the sense of Definition \2. l\ 
Moreover, there exist 5, L > and a semigroup S: x D ^ D such that 

1. V^{u(^u + V-{R]Aq):T\{u-u) <5]. 

2. For all u £ V, Squ = u and for all t,s > 0, StSsU = Sg+tU- 

3. For all u,u £ V and for all t,t > 0, 



\StU — Sf:u\\j 1 < L • i \\u — u\\^i + 



t - t 



4- Ifu gV is piecewise constant, then for t small, Stu is the gluing of Lax 
solutions to standard Riemann problems at the points of jump x ^ 
in u, and of solutions to l^l.3\) - fL^) in the sense of Definition \2.1\ at 
the junction at x = 0. 

5. For all u £ D, the orbit t i-^ SfU is a weak ^-solution to U.l\) . 

6. For any sequences £ Ai, £ Ar and uf € Ui, call the 
corresponding semigroup. If af ^ ai G — > G Ar and 

^ ui £ Ul, then converges uniformly on any compact time 
interval to the semigroup S defined by ai,ar and ui in Lj-^p(M; j4o). 
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Suitable choices of ^ in ()2.3p allow to recover various definitions in the 
current literature. In the examples of the next table, the first component 
$1 of $ is equal to aiqi — a^qr-, ensuring the conservation of mass. 







Meaning 


(L) 


arP{ur) - aiP{ui) 


Conservation of linear momentum, 

see [6j 


(P) 


P{pr)-P{pl) 


Equal pressure, typically motivated 
by static equilibrium, see [21 [3] 


(P) 


P{Ur) - P{ui) 


Equal dynamic pressure, see [Sj [7] 


(S) 


arP{ur) - aiP{ui) 


Limit of the condition for smooth 
variations of the pipes' sections, 
see P[I2] 



The leftmost column displays the letter used below to refer to the solutions 
of (|l.ip - (|1.2p with ^i{ai^ pi, qf, Or, Pr,qr) = drQr—aiQi and ^2 as in the second 
column. The function R in (S) is defined in (12. Sp . see also ()2.1ip . We use 
the above conditions to select particular weak solutions to (jl.ip ~ (|1.2p . 
The main analytical properties of these choices are summarized below. 

Proposition 2.3 All junction conditions <I>: x Aq)"^ — > such that 
^i{ai,pi,qi;ar,Pr,qr) = cirqr-aiqi and ^2 is given by (L), (p), (P) or (S) 
satisfy (^Q), (^^, (^E) and (^E}- Property is satisfied by (p), 

(P) and (S), not by (L). 

Besides, below we characterize the set of initial data ui, Uj. such that ()1.3p 
admits a solution and such that this solution is unique. 

Throughout, we refer to the forward, respectively, backward. Lax curves 
of the i-th family q = Li{p; po, qo), respectively q = L~ {p; po, qo)', see (14. 2p 
for the explicit expressions. As usual, Aj and are the i-th eigenvalues 
and right eigenvectors. The explicit expressions are collected in Section HI 
together with several technical proofs. Moreover, for u (z Aq, the densities 
(pi{u), ^i{u), (pr{u), (Pr{u) denote the /O-component of the intersections of 

the Lax curves p 1-^ (^p, L2 {p;u)^ and p 1— > (yp, Li{p\u)) with the boundary 
of the subsonic region Aq. They are uniquely determined by: 

\2\(pr{u),L:^ {<:pr{u)]u)\ = 0, M Upi{u) , Li {Lpi{u)]u)\ = 0, 
\iyipr{u),L2 {<^r{u)]u)j = 0, Xiyipi{u),Li{(pi{u)-,u)j = 0. 

see Figured! Note that (pi{u) < (pr{u), ^r{u) < ^i{u), (pr{u) < 0riu), 
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p p 

Figure 1: Left, the densities 'pi(u) and <fi{u), along a 1-Lax curve; right, the 
densities <friu) and <fr{u) along a reversed 2-Lax curve; all defined by (|2.4p . 

2.1 (L)-solutions 

Under condition (L), (j2.3p reads 

ai qi = Ur Qr and ai P{puqi) = ar P(pr, Qr) 

and yields the standard Riemann solver giving Lax solutions. This is an 
immediate extension of the standard Lax solutions of (jl.ip . 

Theorem 2.4 Let (EoS) hold and fix positive ai,ar. Consider the Rie- 
mann Problem il.3\) with Ui,Ur £ Aq. Define 

I' = inf |p G [(^;(-u;),(^i(ni)]:aiLi(p;ni) < arL2'((^r(^^r);Wr)| 

= SUpjp G [(^i(u;),(^i(Ui)]:aiLi(/9;Ui) > a,. L2 ('^-r(itr);'"r)| • 

A necessary and sufficient condition for the existence and uniqueness of 
an (L)-solution to attaining values in Aq, is 

aiP{l\Li{l'-ui)) < arP(g(^^Li{r;ui)y^Li{r-ui) 
aiP[l'\Li{l"-ui)) > arP(g(^^Li{r';ui)y^Li{r';ui) 

where g is the inverse in Aq of the map p Ll^ip] u^). 

Note that /' < since the sets defining I' and /" are not disjoint. 

An example of non existence of solutions to (11.30 is provided by the 
following corollary. 

Corollary 2.5 Under the assumptions above, in each of the two cases 

ai < Or and (pi{ui) < (fir{Ur) 
ai > Or and (pi{ui) > (Pr{Ur) 

see Figure\^ an (L)-solution to I attaining values in Aq, does not exist. 



(1 K\ 
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Prom the physical pomt of view, when ai = condition (L) is the 
most reasonable one, since it states the conservation of mass and linear 
momentum. In the case of an elbow, an analog of this condition can be 
justified through the conservation of the linear momentum along a direction 
dependent on the geometry of the elbow, see [6l Propositions 3.2 and 4.2]. 
Por a study of the dynamic of a fluid in a kink, see also [15]. However, in 
the present case of a junction between coUinear pipes with different sections, 
condition (L) is hardly acceptable, see |6t Pigure 7]. 

2.2 (p)-solutions 

Under condition (p), (|2.3|) reads 

ai qi = ar Qr and p{pi) = p{pr) 

and was considered in [H Paragraph 5.2], see also [2j. It was there introduced 
neglecting the pressure drop at the junction, on the basis of engineering lit- 
erature on the subject, see for instance [201 Section 6.3.1]. In this condition, 
the role of the fluid speed is limited to ensure the conservation of mass. 

Theorem 2.6 Let (EoS) hold and fix positive ai,ar. Consider the Rie- 
mann Problem with ui,Ur € ^o- Define 

I' = max \^(p I (ui), ipr{ur)} and I" = mm\^ipi{ui),ipr{ur)} ■ 

A necessary and sufficient condition for the existence and uniqueness of 
a (p)-solution to lll.cl]) attaining values in Aq is 

aiLi{l';ui) > Or L^il'-jUr), . 
aiLi{l";ui) < arL2{l";Ur)- 

An example of non existence of (p)-solutions to (jl.3p is provided by the 
following corollary. 

Corollary 2.7 Under the assumptions above, in each of the two cases 

(pi{ui) > (pr{Ur) or ipi{ui) < (pr{Ur) 

a (p)-solution to U.3\) attaining values in Aq does not exist. 

An example of lack of continuous dependence for (p)-solutions is in [6]. 

^Prom the physical point of view, condition (p) is fully justified in the 
static situation iji = ijr = 0. It is sometimes extended to non static cases in 
the engineering literature, possibly corrected through suitable pressure loss 
coefficients, see [201 Section 6.3.2]. 
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p p 

Figure 2: The situations of Corollaries [131 and 

2.3 (P)-solutions 

Under condition (P), (12. 3p reads 

ai qi = ttr Qr and P{phqi) = P{pr, Qr) 
and was considered in [H [7] . 

Theorem 2.8 Let (EoS) hold and fix positive ai,ar. Consider the Rie- 
mann Problem hl.S^) with Ui,Ur G Aq. Define 

I' = inf |p G i^piiui), i^i{ui)]:ai Li{p; ui) < a,- L^iifriur); 

I" = supjp G [(^i(ni),(^i(ni)]:a;Li(p;ni) > a^Lg ((^r(^ir);^ir)} 

A necessary and sufficient condition for the existence and uniqueness of 
a (P)-solution to i f attaining values in Aq, is 

P{l\L,{l'-,ui))-p(g(^Li{l'-,ui)),^Li{l'-ui)) < 0, 

p(/",Li(r;uO)-^(^5(^^ia";^/)),^^iG";^/)J > o, 

where g is the inverse in Aq of the function p Ll^ip] Ur). 

The proof is similar to that of Theorem 12. 4| hence we omit it. An ex- 
ample of non existence of (P)-solutions to (jl.3p is provided by the following 
corollary. 

Corollary 2.9 Under the assumptions above, in each of the two cases 

(piiui) > (priUr) or (pi{ui) < (pr{Ur), 

see Figure\^ a (P)-solution to !il.'J\) attaining values in Aq does not exist. 

^From the physical point of view, condition (P) implies the conservation 
of linear momentum along directions orthogonal to the pipes, as shown in [71 
Lemma 2.2]. 
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2.4 (S)-solutions 

Under condition (S), (12. 3p reads 

l-CLr 

aiqi=arqr and arP{pr,qr) = aiP{pi,qi) + p {R{a; pi, qi)) da, 

Jai 

where (i?(a; pi,qi), Q{a; Pi,qi)) is the solution to the Cauchy Problem 



£Q{a) = -Q{a)/a, 



3^ aP{R{a),Q{a))\ = p {R{a)) 

Q{ai) = m, 

R{ai) = Pi- 



(2i 



This kind of solution was considered in [9j as limit of solutions to the p- 
system with smoothly varying section. This argument is based on the next 
lemma. 

Lemma 2.10 Fix x < x E M, > and > with ai ^ a,.. Let 
a G C°'i(M;R+) satisfy 



a(x) = a; if X < X 

a strictly monotone if x £ [x, x] 
a{x) = ar if X > x 



(2.9) 



Call R°'{x; pi,qi) the p component of the solution to the Cauchy Problem 



dx{a{x)q) = p{x) 
dx(aix) {q^/p + p{p)fj =p{p)dxa q{x) 



Pi 
m ■ 



Then, the function 

{ai,ar; pi,qi;x,x;a) 



p {R"-{x;pi,qi)) a'{x)dx. 



(2.10) 



is well defined. Moreover, if a}, a? are monotone functions satisfying 112. 9\) . 
then the corresponding functions ^2.10\) coincide. Hence, the following map 
is well defined: 



{ai,ar;pi,pr) 



p{R{a;pi,qi)) da. 



(2.11) 



0.1 
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For a proof of this result, see [9l Proposition 2.7]. Note that the Cauchy 
Problem (j2.8p can be rewritten in the form 

' Q{a) = ^qi, 

da'^y"') a a2p'[R(a))R2(a)-afqf 

R{ai) = pi. 

V 

Lemma 2.11 Consider the system i2.12\) and a point {pi,qi) € Aq. 

1. If qi > 0, then Q{a) is strictly decreasing. If qi = 0, then Q{a) is 
constantly equal to 0. If qi < 0, then Q{a) is strictly increasing. 

2. If qi = 0, then R{a) is constantly equal to pi. If qi ^ 0, then R{a) exists 
and is strictly increasing for every a > ai. Moreover it is bounded for 
a > ai; see Figure\M 




Figure 3: Solutions to (I2.12p starting from various initial conditions. 



Theorem 2.12 Let (EoS) hold and fix positive ai < a^. Consider the 
Riemann Problem with ui,Ur S ^o- Define 

I' = R(^ar;ipi{ui),Li {ipi{ui);ui)^ , 
I" = R(^ar;ipi{ui),Li {(pi{ui);ui)^ , 

where R is the function defined in ^2.12\) . A necessary and sufficient con- 
dition for existence and uniqueness of an (S)-solution to {H^, attaining 
values in Aq, is 



aiLi{0i{ui);ui) > a,- {l';Ur) , 
ai Li {(pi{ui); ui) < Ur {I"; Ur) . 



(2.13) 
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Corollary 2.13 Consider the Riemann Problem 1^1. 3\) with ai < and 
such that (pi{ui) > (pr{ur). Then an (S)-solution to hl.S^) . attaining values 
in Aq, does not exist. 

^From the physical point of view, condition (S) is justified as the hmit 
of smooth changes in the pipes' section, see [l2l Theorem 2]. Viceversa, 
consider n consecutive junctions sited at, say, = i/n, with i = 0, 1, . . . , n, 
separating pipes with section, say, af = a; + {i/n){ar — a/). Imposing condi- 
tion (S) on each junction is equivalent, in the limit n +oo, to the usual 
model for pipes with a smoothly varying section, see |19|, Section 8.1]. 



3 Numerical Examples 

The paragraph is devoted to show results of numerical integrations of solu- 
tions to (|l.ip - (jl.2p with the different choices of the function $ considered 
above. The integrations below are "exact", in the sense that they amount 
to the solutions of Riemann problems at the junction, which are obtained 
through the (approximate) computation of the intersection between Lax 
curves ()4.2p . Neither time steps, nor meshes, are involved. 
Throughout, we fix 7 = 1.4 and the 7-law p{p) = p^. 

(L)-solution - p (P)-solution - p 



(L)-solution - q 



(P)-solution - 



(p)-solution - p 



(p)-solution - q 



(S)-solution - p 



(S)-solution - q 



Figure 4: Here, ai — ar — 2, pi = pr = 1, Qi = Qr = 0. Note that only 
the (L)-solution is not stationary. 

First, consider an initial datum with the same density and zero speed 
in the two pipes: qi = Qr = and pi = pr. Conditions (P), (p) and (S) 
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yield the stationary solution: no transfer of fluid between the pipes and no 
pressure difference at the junction, see Figure [H The (L)-solution signifi- 
cantly differs: it is not stationary, prescribes a transfer of fluid from the right 
(larger) tube to the left (smaller) tube and yields a pressure difference at the 
junction, see Figure [H Indeed, the (L)-solution consists in a 2-rarefaction 
moving to the right and a 1-shock to the left, with fluid flowing from the 
right to the left. 

A situation frequently considered in the engineering literature, see for 
instance [20, Section 6.1.3] and [1], is that of a shock wave hitting a junction. 
We consider the case of stationary flow of fluid moving from the smaller 



{P)-solution - q 



(L)-solution - q 



(p)-solution - p 



(p)-solution - q 



(S)-solution - p 



(S)-solution - q x 



Figure 5: Here, a; = 1, = 2, and the left state is pi = 2.2, qi = 1.387882. For 
the other initial data see (|3.ip . 

pipe, say, on the left, to the larger one on the right. Differently from the 
previous case, we now want to perturb a stationary but non-static situation. 
Therefore, we are bound to choose different stationary configurations for 
the different solutions. Let a; = 1, a,. = 2 and fix the unperturbed state 
pi = 2, qi = 1 on the left. Note that qi > 0, so that fluid flows rightwards. 
Correspondingly, we find the state {pr,qr) to be assigned to the right pipe, 
so that we have stationary solutions in the different cases: 

Solution pr Qr Solution pr qr 

(L) 1.2520452 0.5 (p) 2.0000000 0.5 (3.1) 

(P) 2.2051669 0.5 (S) 2.1064869 0.5 

Obviously, qr is uniquely determined by qi and the sections a^, through 
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mass conservation. Then, a 2-shock moving in the left tube towards the 
junction has left state pi = 2.2, qi = 1.387882 and right state {puQi)- These 
latter values, as well as the forthcoming solutions, are determined using ()4.2p 
and the definitions of solutions. 

The qualitative behavior of the solution to the Riemann Problem (II. 3p 
is the same in all cases: the shock interacts with the junction leading to the 
formation of two waves. A refracted shock proceeding in the right tube and 
a reflected rarefaction moving leftwards, see Figure [5l 

Consider now a stationary configuration perturbed by a shock coming 
from the right (larger) pipe. In cases (p), (P) and (S) we use the same 



(L)-solution - p 



(L)-solution - q 



(P)-soluti( 



(P)-solution - q 



(p)-solution - p 



(S)-solution - p 



(p)-solution - q 



(S)-soIutioii - q 



Figure 6: Here, a; — I, ar = 2. In case (L) pi = 2, qi = 1, pr = 1.2520452, 
Qr ~ 0.5. In cases (p), (P) and (S) pi = Pr ^ 1, Qi = Qr = 0. The interaction of 
a 1-shock from the right results in a 1 shock being refracted in the left tube and a 
2-shock reflected in the right tube. 

data as in Figured! namely ai = 1, ar = 2, pi = pr = 1, qi = Qr = 0, so 
that the fluid is at rest. In case (L), the data a; = 1, = 2, pi = 2, qi = 1, 
Pr = 1.2520452 and Qr = 0.5 yield a stationary solution. We perturb these 
stationary solutions with a 1-shock heading towards the junction from the 
right tube, obtaining the Riemann problem with right data pr = 1.3 and 
(jr = —0.4160525 in cases (p), (P) and (S), and with right data pr = 1-55 
and (jr = 0.1994594 in case (L). Then, the result of this interaction is in 
Figure El In all cases, the interaction results in two shocks, one reflected 
and one refracted. 
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4 Technical Details 



Recall the expressions of the eigenvalues and eigenvectors of the p-system 



ri{u) 



-1 



r2{u) 



1 



(4.1) 



(4.2) 



The Lax curves and their reversed counterparts have the expressions 
Li{p;Uo) = 

L2{p\Uo) = 

Li{p;Uo) = 



Proof of Proposition [K^ A direct application of the Implicit Function The- 
orem allows to define neighborhoods Ai,Ar of ai,ar and Ui,Ut of ui,Ur such 
that for all a; G £ Ar and ui GUi, there exists a unique Uj- G Ur such 

that $(a/, It/; ttr, «r) = if and only if Ur = T{ui; ai,ar). 

The proof of 1.-6. is extremely similar to various results already appeared 
in the literature. More precisely, 1.-5. are shown as in [TJ Theorem 3.3], 
[SI Theorem 3.2], Theorem 2.3]. The latter estimate follows from [8l 
formula (3.2)] and (*l0D. □ 

Proof of Proposition [27^ Concerning ($[0]), in the different cases, simple 
computations show that 



Po 




dr 






'^i P< Po 


po 




- Po) 


{pip) 


-p{po)) 


'^i P> Po 


po 




- Po) 


{pip) 


-p{po)) 


'^i P< Po 


Po 




dr 






P> Po 


po 




- Po) 


{pip) 


-p{po)) 


P< Po 


^ 

po 


20-Pir-^ 


dr 






P> Po 


fo 


lo+^jf-sp- 


- Po) 


{pip) 


-p{po)) 


if P > Po, 


fo 


lo+Pir-? 


dr 






if P < Po, 



Du^^iai,ui;ar,Ur 



—ttr^ \\{Ur) \2iUr) 
ar p'iPr) 

— Qr Xi{Ur) X2iUr) 

— Or Xi{Ur) X2iUr) 



in case (L) 
in case (p) 
in case (P) 
in case (S) 



which does not vanish in Aq. 

In case (S), consider (^[T]). The condition ^{a, pi,qi;a, pr,qr) = im- 
plies qi = qr and so pi = pr] hence (^[T]) holds. Concerning ($[5]) and 
assume that iai, pi, qi; ar, Pr,qr) = 0, which is equivalent to 

I aiqi = arqr, 

I pr = Riar;pi,qi). 
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Since aQ{a) = aiqi for every a and the equation for R depends on afqf, then 
we deduce that, if pr = R{ar; pi, qi), then pi = R{ai; pr, qr) = R{ai; pr, —qr)- 
Thus ^{ai, pi,qi;ar, Pr,qr) = is equivalent to ^{ar, Pr, —qr', cli, pi, —qi) = 
and ($[2]) holds. The fact that aQ{a) = aiqi for every a and the equation 
for R depends on afqf implies that (^(31) also holds. Property ($|4]) follows 
from 2. in Lemma 12.111 

The other cases are immediate. □ 



Proof of Theorem \2.4\ Assume that there exists an (L)-solution u attaining 
values in Aq. Then 

{p{t,0-),Li{p{t,0-y,ui)) G ^0, {p{t,0+),L^{p{t,0+);Ur)) G 

and /' < 0— ), /" > p{t,0—). Moreover, by ^ Lemma 1], we have 

aiP{l',Li{l';ui)) -arP (^g (^^Li(r; n^)) , u/)^ 
< aiP {u{t, 0-)) - QrP {u{t, 0+)) = , and 

aiP (/", ui)) - arP (^g (^^Li{l"; n,)) , ^Li^; ui)^ 

> aiP {u{t, 0-)) - arP {u{t, 0+)) = 

Hence ()2.5p is a necessary condition. 

Assume now that (j2.5p holds. Since I' < I", then 

{p, Li {p; ui)) eAo, ^g (^-^^i {p; n;)^ , ^Li {p; ui)j G Aq 

for every p G [I', I"]- Moreover the function 

p 1 — > aiP {p,Li{p;ui)) - OrP (^g (^^Li{p;ui)^ , ^Li{p;ui)^ , 

defined in [I', I"], is increasing, by P, Lemma 1], and continuous. Hence 
there exists p € [I', I"] such that 

aiP (p, Li{p; ui)) - QrP (^Li{p; ui)^ , -^Li{p; ui)^ = 0. 

Therefore the traces of an (L)-solution at J are given by 

{p,Li{p;ui)) and { g I —Li{p;ui)] , —Li{p;ui) ] . 
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This concludes the existence proof. 

Assume that the Riemann Problem (jl.3p admits two different (L)-solu- 
tions, attaining values in Aq, denoted with {p',q') and {p",q"). Note that 

(?'(t,0-) = Li{p'{t,0-y,ui), q"it,0-) = Li{p"{t,0-);ui), 
(?'(t,0+) = L2-(p'(t,0+);n,), g"(t,0+) = 0+); n,) . 

Without loss of generality, we may assume that p'{t,0+) < p"(t,0+). By [51 
Lemma 1], we deduce that q'{t,0+) < q"{t,0+) and 

arP{p'{t,0+),q'it,0+)) < arPip"{t,0+),q"{t,0+)). 

By definition of (L)-solution, the previous inequality becomes 

aiP{p'{t,0-),q'{t,0-)) < aiP{p"{t,0-),q"{t,0-)), 

and, by [5j, Lemma 1], we deduce that p'{t, 0—) < p"{t, 0—) and so q'{t, 0—) > 
q"{t, 0-), which is in contradiction with q'{t, 0+) < q"{t, 0+). □ 

Proof of Theorem \2.b\ Assume that a (p)-solution {p,q) in exists. Then 

{p{t,0-),Li{p{t,0-y,ui)) G ^0, {p{t,0+),L^{p{t,0+);Ur)) G 
and so I' < 0— ) = 0+) < Moreover, by [5l Lemma 1], we have 

aiLi{l';ui) - Or 
> aiLi{p{t,0-);ui) - Ur ip{t,0+);Ur) = 0. 

aiLi{l";ui) - Or L^if^Ur) 
< ai Li{p{t, 0-); ui) - ar-L^ipit, 0+); u^) = . 

Hence ()2.6p is a necessary condition. 

Assume now that (j2.6p holds. First, we claim that < Suppose, 
by contradiction, that I' > I" . We have two different possibilities: either 
'Piiui) < (friur) or (pr{ur) < <f>i{ui). Consider only the first case, the second 
one being similar. We easily deduce that /' = Lpr-{ur) and /" = ipi{ui) and so 

Li{l';ui) < and Li{l";ui) > (/";'Ur). (4.3) 

Moreover, ()2.6p implies that 

Li{l";ui) ~ Ur ~ Li{l';ui) ' 
which is in contradiction with (j4.3p . Since /' < then 
{p,Li{p;ui)) G Aq, {p,L2{p;Ur)) G 
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for every p € [I' ,1"]. Moreover the function 

p I — > aiLi{p;ui) - ar L2{p;Ur), 

defined in [I', I"], is decreasing, by [SJ Lemma 1], and continuous. Hence 
there exists p £ [I', I"] such that 

ai Li{p; ui) - a.,. L^ip; Ur) = 0. 

Therefore the traces of a (p)-solution at J are given by 

{p,Li{p;ui)) and {p,L^{p;Ur)). 

Assume now that the Riemann Problem ()1.3p admits two different (p)- 
sohitions {p',q') and {p",q") attaining values in Aq. Thus, we deduce that 

p'it, 0+) = p'{t, 0-) / p"{t, 0-) = p"{t, 0+) 
for a.e. t > 0. Without loss of generality, we suppose that 

p'{t, 0+) = p'{t, 0-) < p"{t, 0-) = p"{t, 0+) 
for a.e. t > and so, by using ^ Lemma 1], 

a^Lg (p'(t,0+);nr.) = ai Li{p'{t,0-);ui) > ai Li{p" {t,0-);ui) 

= Or ^2 0+); Mr) > ttj. {p' {t, 0+); Ur) , 

for a.e. t > 0, which is a contradiction. This concludes the proof. □ 

Proof of Corollary \2.7\ Consider only the case 0r{ur) < ifi{ui), the other 
one being similar. Assume by contradiction that {p, q) is a (p)-solution 
to (II. 3p attaining values in Aq. Therefore, we deduce that 

q{t,0-) = Li{p{t,0-);ui), 
q{t,0+) = (p(t,0+);ur.), fora.e. t>0, 
p{p{t,0-)) = p{p{t,0+)), 

hence p{t,0—) = p(t,0+) for a.e. t > 0. Moreover, p{t,0—) > (pi{ui) and 
p{t, 0+) < (pr{ur), which gives a contradiction and proves non existence. □ 

Proof of Corollary \2.9l Consider only the case (pr{ur) < (pi{ui), the other 
one being similar. Assume by contradiction that (/j, q) is a (P)-solution 
to (jl.3p . attaining values in Aq. Therefore we deduce that 

(?(t,0-) = Li{p{t,Q-)-ui), 
(?(t,0+) = L2(p(t,0+);M,), 
P(p(t,0-),g(t,0-)) = P(p(t,0+),g(t,0+)), 
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for a.e. t > 0. Moreover we have 



P{p{t,0-),q{t,0-)) > P{ipi{pi,0),Li{ipi{pi,0);ui)) 

> Pip{t,0+),q{t,0+)) 
and so we get a contradiction, proving not existence. □ 

Proof of Lemma \2.11[ The first item is trivial; so we consider only the 
second one. 

If qi = 0, then ^R{a) = and R{a) is constantly equal to pi. Assume 
therefore qi ^ 0. The sign of the derivative of R{a) is given by the sign of 

ay{R{a))R\a)-aUl 

At the initial point a;, this term is {p'{pi)p'i ~ 9^)' '^^ich is strictly pos- 
itive if and only if {pi,qi) € ^o- Thus R{a) is increasing in a neighborhood 
of ai. If R{a) is increasing, then we deduce that 

a'p'iR{a))R\a) - afqf > af\p'{pi)pf - qf 

for every a > ai, provided R{a) exists. Using this estimate and the com- 
parison theorem for ODE, we conclude that R(a) exists and is increasing 
for every a > ai. Finally, for a > ai sufficiently big, the derivative R'{a) 
can be bounded by K/a^, where K is a constant depending on the initial 
conditions. Thus we deduce that R(a) is bounded. □ 



Lemma 4.1 Fix {pi,qi) € Aq and denote with R{a;pi,qi) the solution to 
the ODE in ^KTE) . Define 



d d 
zi{a) := —R{a;p,qi)\p=p^ and Z2{a) := —R{a] pi,q)\g=g^. 

Then z\ and satisfy the following system 



zi{a) = {a,R{a;pi,qi),qi)) zi{a) 



A", — 9 

da 

^Z2{a) = ^g{a,R{a;pi,qi),qi)) Z2{a) + -§^g{a,R{a;pi,qi),qi)) 
zi{ai) = 1 
Z2{ai) = 0, 

, . ^ P afq^ 
where g (a, p,q) = — 



a a^p'{p)p'^ — afq"^ 
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The proof consists in the classical derivation of the solution to an ODE 
with respect to a parameter, hence we omit it. 

Proof of Theorem \2.12l Consider the following curve on 
■0: [01(^1,91), 'fi{pi,Qi)] — > Ao 

(4.5) 

where ^pl{s) = R{ar; s, Li{s; pi,qi)) and ^2(5) = ^Li{s; pi,qi)). Since the 
point (s, Li(s; pi,qi)) belongs to Aq for every s in the domain of ip, then 
we easily deduce that the image of tp is contained in by Lemma [2.111 
Clearly the second component of ■0 is decreasing with respect to s. 

We claim that ipi is increasing with respect to s. In the following we use 
the same notation of Lemma 14.11 The derivative of ipi is 

i^'iis) = ^ R{ar;s,Li{s;pi,qi)) 
as 

d _ _ d _ _ d 

= -Q^R{ar;s,Li{s;pi,qi)) + — R{ar; s, Li{s; pi, qi))— Li{s; pi, qi) 

_ _ d 

= zi{ar; s, Li{s; pi,qi)) + Z2{ar; s, Li{s; pi,qi))-^ Li{s; pi, qi). 

By ()4.4p . we have that zi{ar; s, Li{s] pi, qi)) > for every s in the domain 
of the curve ip. Let us consider some different cases. 

1. Li{s, pi, qi) > 0. In this case the derivative 

d 

g^aiar, R{ar;s, Li{s, pi, qi)),Li{s, pi,qi)) 

is strictly positive. For a > ai, define the function 

_ _ d 

P(a) = zi{a;s,Li{s;pi,qi)) + Z2{a; s, Li{s; pi, qi))— Li{s; pi, qi). 
Easy computations show that 

= §-p9{a,R{a]S,Li{s,puqi)),Li{s,puqi))(3{a) 
+ -§^9{a, R{a; s, Li{s, pi,qi)), Li{s, pi, qi)), 
P{ai) = 1. 

Define 

a = inf {a > a/ : /?(a) = O} . 

Assume by contradiction that a < +00. In this case we deduce that 

d _ d _ _ 

= -g^g{a,R{a;s,Li{s,pi,qi)),Li{s,pi,qi)) > 

by assumptions and this is not possible. Hence /3(a) > for every 
a > ai. In particular P{ar) > and so V'i(^) > 0- 
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2. Li{s, pi,qi) = 0. In this case the derivative 

d 

—g{ar,R{ar;s, Li{s, pi, qi)),Li{s, pi,qi)) 

vanishes and so, by ()4.4|) . Z2{ar; s, Li{s; pi, qi)) = 0. Hence ip'i{s) > 0. 

3. Li{s, pi,qi) < 0. In this case the derivative 

d 

— c/(ar., Riar] s, Li{s, pi, qi)),Li{s, pi,qi)) 
is strictly negative and so, by (j4.4p . Z2{ar; s, Li{s; pi, Iji)) < 0. Hence 

> 0. 

By the previous considerations, a necessary and sufficient condition for the 
existence and uniqueness of an (S)-solution {p, q), attaining values in Aq, is 
that the image of ^ intersects in a unique point the image of the curve 

['fr{Ur),'fr{Ur)] > Aq 

( T-i - <\ 4.6 



The image of the curve (14. 6p divides the set in two parts and this permits 
to conclude. □ 

Proof of Corollarv \2.1'A Assume by contradiction that (p, x) is an (S)- 
solution to (II. 3p . attaining values in ^o- Therefore we deduce that 

g(t,0-) = Li(p(t,0-);720, 
g(t,0+) = L2(p(t,0+);n,), 
ii(a,;p(t,0-),Q(t,0-)) = p(t,0+), 

for a.e. t > 0. By Lemma [2.111 we deduce that i?(ar; /o(t, 0— ), 0— )) > 
0— ). Moreover, by hypotheses, we have that 

P(t,0-) > (pi{ui) > ipriPr,qr) > p(t,0+) 

and so we obtain a contradiction. □ 
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